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ABSTRACT 


New long path,, low temperature , moderate resolution spectra of 
methane and ammo nia, broadened by hydrogen and helium, are used to cal- 
culate non-gray model atmospheres for the four jovian planets. Solar 
energy deposition in the upper atmospheres is of major importance; the 
fundamental and' first overtone of hydrogen contributes enough absorption 
to create a thermal inversion for. each of the planets . A wide range of 
models of emergent flux, dependent on assumed mixing ratios and scattering 
by condensates and photochemical organics, is used to calculate the suite 
of emergent spectral fluxes and representative limb darlcenings and 
brightenings in the range 600 to 12,000 cm" ^ for comparison with the 
Voyager infrared spectra. The temperature differences between jovian 
belts and zones corresponds to a difference in the ammonia cirrus particle 
radii (1 to 3 ym in. zones; 10 ym in belts). The jovian tropopause is ap- 
proximately at the 0.1 bar level. A thin ammonia cirrus haze should be 
distributed throughout the Saturnian troposphere; and NHg gas must be 
slightly supersaturated or ammonia ice particles carried upwards convectively 
in the upper troposphere of Saturn. Substantial methane clouds exist on 
both Uranus and Neptune, implying that considerable caution should be used 
in the application of simple reflecting models to interpret the spectroscopy 
of both planets. There is some evidence for almost isothermal structures- 
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in the deep atmospheres of Uranus and Neptune. Expected seasonal 
variations during the Uranian year should be detectable in the emergent 
flux. There is either an internal heat source on Neptune or an abundance 
of particles with radii larger than 10 pm. 



Silvaggio (1977) has described a new set of observations of 
the absorption of methane and arnronia in the infrared, using a multiple 
reflection folded path 'White cell with path lengths up to 72 m. Tem- 
peratures were varied between about 110° K and about 275° K and the 
effects of foreign gas band broadening by hydrogen and helium were 
studied over this temperature range. A subset of these observations 
at a resolution between 0.25 and 1.0 cnT^ measured over the spectral 
range 600 to 12,000 cm - '** lend themselves to the construction of model 
atmospheres and emergent fluxes for the four jovian planets, Jupiter, 

Saturn, Uranus and Neptune. 

hi order to reduce the number of frequency intervals needed for 
structure calculations , we used average monochromatic opacities 
(cf. Pollack, 1969). Because of atmospheric absorption in the visible 
and near infrared, each atmospheric layer has a non-zero absorption of 
solar energy and there will be a variation of the net flux with altitude. 
Effects due to the general circulation of the atmosphere are neglected. 

Our non-gray atmospheric models are based on the following additional as- 
sumptions : > He, GH^, MLj and possibly EpO are the chief sources of 

Infrared opacity. Scattering is approximated either by Rayleigh scattering 
or by Mie scattering for spherical particles . The atmosphere is in radia- 
tive equilibrium, except when the temperature gradient exceeds the adiabatic 
gradient, in which case the temperature gradient is set equal to the 
adiabatic gradient. The atmosphere is in hydrostatic equilibrium. The 
unsaturated components are conpletely mixed. The plane parallel atmosphere 


- 1 - 



- 2 - 


approximation holds and the acceleration due to gravity is constant with 
altitude. The populations of the molecular energy levels are given by 
local thermodynamic equilibrium. The heat capacity of the atmosphere is 
large enough and the planetary rotation is fast enough that the effective 
temperature at a given latitude is independent of longitude. There is an 
opaque "surface," probably a sharp cloud deck, -with a frequency-independent 
reflectivity. Saturated components , 1^0, and CH^, are assumed to form 

clouds, and the quantity of condensates is derived from the difference be- 
tween the saturation vapor pressure and the partial pressure implied by 
the equation of hydrostatic equilibrium. 

We solve the equation of radiative transfer iteratively. With an 
initial temperature profile and gas abundances for each constituent, a 
model atmosphere is generated, allowing for convection. The atmosphere is 
quantized, creating a finite number of layers and, for the purposes of 
radiation balance, a number of angles whose values are determined as those 
for Gaussian quadrature. For the conditions present in each layer, absorp- 
tion and scattering coefficients are calculated and the equation of radiative 
transfer solved. The new temperatures derived in this manner replace the 
previous estimate, and the procedure then iterates until the temperature 
profile undergoes an arbitrarily small change. The emergent flux is con- 
strained to equal the observed value. The temperature structure of the at- 
mosphere is approximated by a series of layers, each of constant temperature 
T, for which the ideal gas law and the equation of hydrostatic equilibrium 
apply. 
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Some atmospheric constituents saturate as the temperature decreases; 
their densities can be approximated as n sat = (n 0 /T) exp [A - B/ (T - C) ] , 
where the parameters- have the values- of Table- 1. 


Table 1 

Vapor Pressure Temperature Reference 

Constants Range 


Gas 

n 

-3 ° -1 

[cm J (K ) '■*•] 
0 

A 

B 

C 



nh 3 

4.118 x 10 17 

25.88 

3752.6 

- 

all 

Lasker (1963) 

CH. 

4 

9.673 x 10 18 
21 

7.35 x 10 

16.5629 

15.4172 

939.94 

933.61 

4.66 

5.38 

T< 60 K 
60 K < T< 90K 
T > 90. K 

Wilholt (1971) 

h 2 ° 

21 

7.35 x 10 

13.9787 

17.45915 

5202. 

6159. 


T > 273 K 
T< 273 K 

Handbook of Chem- 
istry and Physics 


When the tsxperature gradient exceeds the adiabatic gradient, the tempera- 
ture gradient is set equal to the adiabatic gradient: 


dT 

B2- 


- -s/c . 


ad 


C P = i n i c P ,i /ai i 


where C . is the specific heat at constant pressure for gas i, and n. is 

P , i 

the density of gas i. Errors in will be small (Dubisch, 1974) compared 

to errors in the ideal gas law. C . is calculated for each molecule by 

P 

C = c, + c 9 T -f -CoT^ cal/deg-mole, 


where the parameters have the values given in Table 2. 
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Table 2 

Constants for Specific Heat 


C-as c x ° 2 

H, 4.387 1.45 x 10' 


He 

4.97 



-2 

CH, 

4 

3.-38 

1.79 

X 

10. 

-3 

nh 3 

5.92 

8.96 

X 

10 

-3 

h 2 ° 

6.89 

3.28 

X 

10 



Reference 
Gopal (1966) 

Zemansky (1957) 


The value for is the specific heat for an equilibrium mixture (Trafton, 1965) . 

Where there is condensation, the adiabatic gradient can be derived 
from the Clausius -Clapyeron Equation, rewritten as 



where H is the heat of vaporization/siiblimation, the average specific 
heat for the unsaturated gases, n the condensate number density, n the 
atmospheric number density, m the mean mass of atmospheric molecule, and m c 
the mass of the condensing molecule. 

The equation of radiative transfer is 


dl /dr 
v v 




where I is the specific intensity, J v the mean intensity, and the 

optical depth. Multiplying both sides by exp (-x ) , integrating over x y 

and assuming J constant over an interval At , we find 
a v \) 


. , t -i -2At -2At 

I 1+1 » I 3 ;" 1 e v + (1 - e v ) J' x 

V V x ' V 


.for the -downward intensity- at level i + 1, in terms* of quantities at 
levels i and i - 1, each of optical thickness Ax^ . A similar relation 
applies for the upward intensity, I^~ i- s calculated by inserting 

the I 1 found in the previous step into the expression 


rr w t 

- JJtt + -~e v P(fl)I Q + (1 - 5) 0 ) B 




o' v 


Oiere m 0 is the single scattering albedo, B y Is the Planck function, 
and I Q represents incident intensity. The phase function P(Q) is con- 


strained so that 


JJ P(J2)dft = 4 tt. 


In the initial computation, the intensity is assumed to be 
zero in obtaining J. The downward intensity is first calculated be- 
ginning at the top of the atmosphere, and the calculations progress 
downward. Then, at reaching the bottom of the atmosphere, the upward 
intensity is calculated. This process is repeated, using the most 
recent estimates of I x , until I converges to a constant or slowly 
fluctuating value. The middle term in Eq. (2) is the source function 
due to the first scattering of the incident radiation. If this term 
becomes comparable to the others, then higher orders of scattering must 
be taken into account by adding the source function for higher orders 
of scattering 

J n + 1 = 2GF p (^) I n dfie " T for (n + 1)* order. 


(3) 
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* 1 

The working equation (2) and its analog for I^ -i can be solved 
interatively. By adjusting the temperature of each level so that 
the radiation flux absorbed equals the flux reemitted, a tanperature 
profile is derived. The correction to the temperature is AT ~ (A-E)/4oT^ 
where A is the flux absorbed, E the flux emitted, and a the Stefan- 
Boltzmann constant. With an initial temperature guess for each layer, 
the program calculates a model atmosphere, then solves the equation of 
radiative transfer and makes a correction to the initial temperature 
estimate. On successive iterations the temperature corrections diminish 
until a solution is reached. We assume that only solar radiation is 
incident from above, and that the lower cloud surface is a partially 
reflecting layer at which some fraction of the incident energy is ab- 
sorbed and the remainder reflected. 

Trafton (1967) and Pollack and Ohring (1973) have previously 
calculated model atmospheres for Jupiter. We will employ our computa- 
tional procedures to obtain improved atmospheric models, in particular 
\ 

by allowing for the temperature- dependence of the ammonia and methane 
opacity. The opacity of the pres sure- induced transitions of hydrogen 
are determined essentially in the same manner as Trafton (1967) . However, 
as Pollack and Ohring (1973) have done, we too allow for the deviation 
of the line profile in the far wings from a Lorentzian profile. Fol- 
lowing a similar technique, the opacity due to the hydrogen fundamental 
vibrational and first overtone band is taken into account. Using the 
same method of calculation as Trafton, it is found that there is good 
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1 

agreement (within 37.) with our observations of H 2 in the 500-1000 cm 
range for, 20 amagats of .H 2 at-- temp eratures of 120K, 195K, and 273K 
(cf.- Fig. I). Theory and observation for the first overtone band are 
compared in Fig. 2. 

In general, ammonia has a broad opacity in the 8 - 14 yin range 
and dominates there, while hydrogen is the chief source of opacity be- 
tween 14 and 40 ym, although there- can be contributions from the pure 

rotation spectrum of water and ammonia. Methane dominates in the 6 - 8 ym 

* 

range, and both methane and ammonia are strong sources of opacity short- 
ward of 3.5 ym. The inclusion of methane and ammonia is important for 
Jupiter , while ammonia is of lesser importance for the other giant planets 
because their .low effective temperatures prevent significant amounts of 
IIH^ from appearing in the upper layers. An entirely empirical expression 
for the self -broadening and foreign gas broadening - as well as temperature 
dependence of methane and ammonia absorption, derived from the White cell 
measurements (Silvaggio, 1977), was used for the model atmospheres. 

We have allowed in our calculations for Rayleigh scattering by 
the molecules in the free atmosphere and MLe scattering by the cloud 
aerosols. We approximate particles of arbitrary shape and composition 
by isotropic homogeneous spheres, because randomly oriented irregularly 
shaped particles should statistically approximate spherical symmetry 
(Sagan and Pollack, 1967) . The spherical Bessel functions of Mie theory 
and their recursion relations were calculated by the method described, by 
Wickramasinghe (1973) . The gas amounts in the layers for our models 
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never exceeded by more than a factor of 2 the range of conditions run 
in the laboratory. Values of the absorption coefficient for methane 
at wavelengths shorter than 0.8 ym, important for calculating the solar 
energy deposition, were obtained from Giver (1978) . 

With a model identical to that in Pollack and Ohring (1973) , 
this technique converged on a temperature profile that was virtually 
the same as their profile. Model . atmospheres having the same properties 
but a different number of levels, frequency intervals or angles have been 
compared. A modest number of frequency intervals will provide an accurate 
model profile. However, the high spectral resolution for the emergent 
flux using the low resolution model profile is only as accurate as the 
low resolution temperature profile. For most cases, a modest number of 
layers, about 20, angles, about 12-16, and frequency intervals, approxi- 
mately 100 cnf ^ in width, give a model that is accurate to within 0.5° K 
for each layer. The improvement by increasing the number of frequency 
intervals' is due to a better averaging of the absorption coefficient over 
the internal, when weighted by the Planck function. 

Tin nature, a distribution. of particle size is usually encountered. 

If n(r)dr is the number of particles per unit volume with a radius be- 
tween r and r + dr, the scattering efficiency becomes 

%ca ' I fr 2 <W r >“<»dr /T r 2 n(r)dr, (4) 

-'O "o 
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where G is the geometric cross-section of the particles. A similar 
expression for C^. applies for the extinction efficiency.- -For -our 
calculations we use the distribution aroloyed by Hansen (1971) , 

n(r) <x r^“^^exp(-r/ab) , (5) 

because it has simple properties and approximates several other distribu- 
tion functions used in the literature, e.g. by Khrgian (1961) and 
Deinnendjian (1964) . Figure 3 illustrates the distribution for several 
values of a and b. 

The values of the index of refraction used for the gases were 

2 

approximated with the Cauchy formula n r -l = a(l + BA ) , A. in ym, with 
coefficients given in Table- 3 . 



TABLE 3 


Cauchy Coefficients 

for Real Part of 

the Refractive Index of Gases 



a 

JU 

H 2 

13.58 

x 10~° 

7.52 x 10“ 7 

He 

3.48 

x 10" 5 

2.3 x 10~ 3 

CH. 

4 

44.1 

x 10~ 5 

- 

nh 3 ' 

37.0 

x 10~ 5 

1.2 x 10~ 3 

h 2 o 

2.5 

-4 

x 10 

- 
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For particles we assumed n r = 1.33 and 1.4, respectively, for H2O and 
CH^. These values do not take any dependence on temperature or wave- 
length into account. The 1^0 clouds are deep in the troposphere of 
Jupiter and Saturn and should have little effect, on energy balance . A 
complete treatment of H^O can be found in Irvine and Pollack (1968) . 

The CH^ particles will be a significant factor on Uranus and Neptune, 
if there is condensation. However, no reliable data have been found for 
CH^ particles in the middle infrared wavelengths . The values for the 
coup lex index of refraction of solid NHg were taken from Taylor (1973) . 

These values were determined from a thin film of NHg at 77°K and do not 
account for any temperature dependence. However, the wavelength dependence 
has been resolved for 10 cnf ~ intervals. For models employing the Sg/organic 
aerosol produced in Jovian simulation experiments by Khare and Sagan (1973) , . 
the values of n(A) were calculated assuming MLe scattering; the results 
should match their published spectrum at room temperature . 

In the present study, a steady-state distribution of temperature 
and density for the atmosphere of a giant planet above an opaque surface 
has been calculated for a variety of models. The variables involved in 
such a study include the brightness temperature, gas abundances, spectral 
dependence of the albedo , and nature and location of Mie scatterers . 

With a very large number of possible combinations in a multi-dimensional 
phase space, and a finite budget, no attempt can be made to explore all 
effects on a fine grid. Father, our approach is to develop a series of 
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models by successively adding an atmospheric constituent and examining 
its effects. We assume that, except for hydrogen, helium, methane, and 
ammonia (and ammonia and methane particles if condensation occurs) , any 
other constituents will have an insignificant effect on the structure 
even though they may be noticeable in the emergent flux of the planet . 

Our program does not require the specification of atmospheric 
temperatures at some lower boundary. Because our only constraint is the 
observed brightness temperature we do not explicitly postulate an in- 
ternal heat source. In the case of Jupiter and Saturn, where the ob- 
served effective temperature is well above the temperature if there were 
no internal source, the temperature of the lower layers were increased 
until the calculated emergent flux equalled that observed. For Uranus 
and Neptune, where there is some doubt as to ’whether there is an internal 
heat source, the bolometric Bond albedo was chosen to be consistent with 
the effective tanperature. The atmospheres of Uranus and Neptune at 
large optical depths may become isothermal below the CH^ cloud (cf . Sagan, 
1969) . This effect can be altered by having larger CH^ ice particles 
effectively scattering the long wavelength radiation. 

Jupiter Atmospheric Structure 

Jupiter is the nearest, brightest, closest and most visited of 
the outer planets, and is accordingly ■ the best understood . Yet, there 
still exist uncertainties in the methane and. airrnonia abundances and there 
is a wide selection of minor constituents from which to choose . The 
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flights of Pioneers 10 and 11 (Carlson and Judge, 1974) and the 0 Scorpii 
occultation results (Elliot et al. , 1974) indicate a helium abundance of ap- 
proximately 10% and an effective temperature of 127° K (Orton, 1975a) . The 
line of sight abundance of CE^ is in doubt, with recent estimates 'ranging from 
50 m-amagat (deEergh et al. , 1976) to 250 tn-amag (Lutz et al. , 1976) . 

Similarly the abundance of NHg ranges from 13 m-amagat (Mason, 1970) to 
30 m-amagat (Woodman et al. , 1977) . The total abundance of seems to be 65 
km-amagat (Margo lis and Hunt, 1972; Fink and Belton, 1969) . It is possible 
that the various amounts of CH^ and NH^, derived by observation of dif- 
ferent bands and assuming a simple reflecting layer, are consistent if 
scattering is treated correctly. 

Sagan and Salpeter (1976) have suggested a wide variety of 
possible particles in the troposphere of Jupiter. A worthwhile treatment 
of scattering requires an estimate of complex refractive indices, size 
distribution and number density. We shall assume that these particles 
will have little effect 'on the energy balance. However, they may be 
noticeable in the emergent flux. For the sake of definitiveness, the 
tropospheric particles will be assumed to be similar to the Sg/ organic 
material synthesized by Khare and Sagan (1973, 1975). In addition, ethane 
is included with the absorption coefficients taken from Khanna et al . 

(1974). The size distribution for the particles [Eq. (5)] was arbitrarily 
but not implausibly chosen to be a = 1, b - 0.01. 
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Results on the Jovian atmospheric structure are shown in 
Figures 4 and 5 for the nodels. of Table 4. The asymptotic temperature 
in the mesosphere is controlled almost exclusively by the abundance of 
methane. ' The minimum temperature at the top of the troposphere is con- 
trolled by the abundance of. anmonia and the size of the NH^-ice particles. 
The larger (10 ym) radii particles scatter the thermal radiation more 
effectively than the smaller (3 or 1 ym) particles. There is no noticeable 
difference between the 1 and 3 ym size particles in the resulting structure. 
No significant difference is found between Models 3 and 7. 

The difference between belts and zones, while not explicitly 
explored, can be appro ximat ed by the large or small particle size models . 

If downward motion effectively removes the NH^ cloud, the small particle 
radii models should accurately represent the cloud- free zones . The. 4° K 
difference between Model 3 and Model 6 is of the same magnitude as the 
3° K variation Qrtcn (1975a) found between the South Tropical Zone and 
the South Equatorial Belt. 

Also shown in Figure 4 are the model atmospheres of Orton (1975a) 
and Sagan and Salpeter (1976) . The low temperature opacities derived in. 
the laboratory raise the tropopause to approximately the 0.1 atm level, 
higher than in the other models. However, the Sagan and Salpeter model 
matches the present results reasonably well. The mesospheric tsnperatures 
in this work do not rise as quickly as those of Orton. 
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Table 4 

Jupiter Model Atmospheres (cf. Figures 4 and 5) 




He 

CH, 

NH_ 

NH 3 Ice 

T 


2 

km-amagat 

km-amagat 

4 

m— amagat 

3 

m-amagat 

r 

avg 
' (pm) 

e 

K 

1 . 

65' 

7.3 

0 

0 

— 

127 

2. 

65' 

7.3 

50 

0 


127 

3. 

65 

7.3 

50 

13 

3 

127 

4. 

65 

7.3 

250 

13 

3 

127 

5. 

65 

7.3 

50 

13 

1 

127 

6. 

65 

7.3 

50 

13 

10 

127 

7. 

65 

7.3 

50 

13 

3 

127 


Other 


C-H-, Khare-Sagan 

L O 

material (2pm radius' 


-15- 


Satum Atmospheric Structure 

While the- theoretical problens of the atmosphere- of. .Saturn, 
are similar to those for Jupiter, our knowledge of Saturn is much -less 
extaisive than that, of Jupiter. • H2 was first identified on Saturn by Munch 
and Spimrad (1963) . The common technique is the observation of the 
4-0 and 3-0 overtones in the quadrupole spectrum. Abundances 
range from 190 km-amagat (Owen, 1969) to 80 km-amagat (Trafton, 1971) . 

Recent estimates of methane abundances range from 35 m-amagat (Trafton, 

1971). to 300 m-amagat (Lutz et al, 1976). The abundance of NHg 011 Saturn 
has been a matter of debate. A theoretical interpretation of the microwave 
spectrum of Saturn by Guilds et al. (1969) requires MHg ; 2 m-amagat are .ob- 
tainedby Woodman et al. (1977). Whether there is enough armonia to be de- 
tected by optical wavelength spectroscopy will depend upon details of the. 
atmospheric structure. There -are no definitive observations for the 
presence of helium. However, it sesns unlikely that if helium is oresent 
its nuudng ratio would differ greatly from that on Jupiter. 

The infrared brightness, temperature of Saturn has been measured 
and we shall use the broad-band measurement of Aumann et al. (1969) of 
97°’ K as a representative value. If the bolometric albedo is 0.61, then 
this effective temperature implies that Saturn, like Jupiter, radiates 
more energy than it receives from the sun, about 3.5 times as- much. 

(See also Erikson et al, 1978) . 

The model of atmosphere of Table 5 is exhibited in Figures 6 and 7. 
Ihe addition of He or MH^ warms' Saturn’s upper atmosphere only slightly, 
while CH^ dramatically raises the temperature of the mesosphere. Because 
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TASLE 

Saturn Model Atmospheres 



H 2 

He 5 

CH. 

4 

km-amagat 

km-amagat 

m-amagat 

1 . 

100 

0 

0 

2. 

200 

0 

0 

3. 

100 

25 

0 

4. 

100 

12.5 

0 

5. 

100 

12.5 

35 

6 . 

100 

12.5 

300 

7. 

100 

12.5 

100 

8. 

100 

12.5 

35 

9. 

200 

25 

300 

10. 

100 

12.5 

300 


5 


(cf. Figures 

6 and 7) 

NH 3 Ice 


NH . 

r 

T 

3 

avg 

e 

m-amagat 

(ym) 

°K 

0 

- 

97 

0 

- 

97 

0 

- 

97 

0 

- 

97 

0 

- 

97 

0 

- 

97 

0 

- 

97 

sat 

3 

97 

sat 

3 

97 


sat 


3 


93 
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of tifae lew temperatures there can be little ^ the -upper atmosphere, 
so it is not --unexpected that -NH^ contributes little to the -absorption 
of sunlight. The amount of GH^ affects the degree of -warming, although 
the sensitivity to mixing ratio is not nearly so large as on Jupiter. 
These results are stunn e r ized in Table 6. 


TABLE 6 

Upper Atmospheric Temperatures on Saturn 
(Maximum/Minimum °K) 


He added 

CH^ (35 m-amagats) 
CH^ (100 m-amagats) 
CH^ (300 m-amagats) 


73/70 

76/73 

112/75; with NH 3 113/75 

115/76 

121/81 


All the models are consistent with a cloud deck in the 120 - 140° K 
range. When is included, a thin NH^-ice haze is distributed throughout 
the troposphere. The reflecting layer could be consistent with the top of 
the thicker portion of an cloud. If the effective tanperature of Saturn 
is as low as 93° K, then the reflecting layer is in the 110° K region. The 
upper atmosphere has a minimum of 70° K and warms to about 108° K. 
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Uranus and Neptune Atmospheric Structure 

Hydrogen has been identified on Uranus and Neptune with the 

observation of the 3-0 and 4-0 bands. The abundances of are 

uncertain because of uncertainty about conditions under which the 

spectral lines are formed and the limited number of observations. 

Estimates range from 200 km-amagat by Poll (1971) to 500 km-amagat by 

Trafton (1974) for Uranus. Quadruoole lines at 1.4^an have' been reported on 

Neptune; the abundance above sens high altitude is quoted as 5 km-a m a g at 

(Joyce et- - al . , ;1977) . CH^ has been observed and the range of abundances in 3 km - 

amagat by Owen (1967) to 15 km-amagat by Cess and Ramanathan (1972) for 

Uranus. Inhomogenous scattering models suggest the low end of this range 

(Macy et~al. , 1978) . While estimates for CH^ on Neptune (Danehy et al. , 1978) 

are 157o mare than on Uranus , the same mixing ratio implies the abundance 

to be a factor of 2 greater for Neptune. (See also Macy et al. , 1978) . 

Ammonia has not been identified in optical spectra of Uranus or Neptune, which is 

expected because of the low temperatures . There is no evidence on the presence or 
absence of He. 

Only a limited number of brightness temperature measurements 
exist for these two outer planets. For the purpose of this paper, we shall 
adopt an equilibrium temperature of 58° K for Uranus and 48° K for Neptune 
(Danielson, 1977) . The large 98° inclination of Uranus must result in 
severe seasonal changes. During 1966, the axis of Uranus was in the plane 
normal to the solar vector, while in 1985 the north pole will face the sun 
continuously. As a result, the effective temperature may then rise to 
63° K assuming there is no effective heat exchange between, hemispheres. 

With these low temperatures , methane clouds should form; for lack of a 
better estimate, the CH /t ice size distribution function parameters will be 

“T 

taken as a = 1, b = 0.01. 



-19- 


The model atmospheres for Uranus and Neptune shown in Table 7 
lead to' the structures exhibited in Figures 8 and 9. Atmospheric 
tsuperatures near the 10 mb level are summarized in Table 8. An es- 
timate of the possible long-term secular dependence of ' atmospheric 
structure on Uranus with orientation of the rotation axis can be gar- 
nered by comparing models 5 and 6 in Figure 8. The addition of helium 
and/or methane to the atmosphere results in a substantial wanning of 
the upper portions of the atmosphere. However, extensive methane clouds 
mast exist on both planets and it is impossible to attribute the large 
CH^ absorption features to the gas above the reflecting layer. Most of 
the observed absorption must be the result of the cloud, either by 
multiple scattering and/or absorption by the methane ice itself. In the 

models considered, the clouds are extremely thick, with column densities 
_o 

> 50 g cm , and it seems very plausible that the troposphere exhibits 
extensive multiple scattering. It is possible that the dramatic shallowing 
of the temperature gradients exhibited deep in the atmospheres of Uranus 
and Neptune represents an approach to isothermal deep atmospheres (in the 
absence of an internal heat source) as described by Sagan (1969). However, 
our calculation techniques are not designed for large optical depths and 
we do not here place extreme confidence in our deep atmosphere results. 
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TABLE 7 


Uranus Models 



h. 

He 

ch 4 

CH^ ice 

T 

e 


km-amagat 

km-amagat 

km-amagat 

r (urn) . 

avg 

°K 

X. 

250 

0 

0 

- 

58 

2. 

500 

0 

0 

- 

58 

3. 

250 

28 

0 

- 

58 

4. 

500 

56 . 

0 

- 

58 

5. 

250 

28 

sat 

3 

58 

6: 

250 

28 

10 

3 

68 


Neptune 

Models 






h 9 

He 

CH 4 

CH 4 ice 

T e 


km-amagat 

km-amagat 

km-amagat 

r avg(ym) 

°K 

1 . 

500 

0 

0 

- 

48 

2. 

500 

56 

0 

- 

48 

3, 

1000 

0 

0 

- 

48 

4. 

1000 

120 

0 

- 

48 

5, 

500 ' 

56 

sat 

3 

48 




TABLE 8 

Upper Atmospheric Temperatures on- Uranus and Neptune 
(Maximum/Minimum °K) 



H 2 

+He 

+CH-. 

4 

Uranus 

30/28 

48/44 

70/48 

Neptune 

2S/23 

30/27 

60/39 


The warm Uranus model appears very much the sane as the normal 
model, except that the troposphere is closer to the adiabatic lapse rate. 
For ice particles substantially larger than 3 inn, the contribution of the 
cloud becomes increasingly important because thermal radi a tion will no 
longer be able to penetrate a thick cloud. 
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Jupiter Emergent Flux 

The modeling technique used for the solution of the radiative 
transfer problem in the atmospheres of the outer planets lends itself 
to the calculation of the brightness temperature as a function of 
wavelength and phase angle. The resolution of the infrared spectrum 
is limited only by the resolution of the absorption and scattering co- 
efficients. Except for Jupiter , the. outer planets have been observed 
primarily over broad band regions of the spectrum. By comparing the 
observed values with our calculated brightness temperatures , the range 
of parameters for the many variables can be narrowed. In this comparison, 
the resolution for the calculated brightness temperatures was kept ap- 
proximately the same as for the observations . Below we propose a set 
of model parameters for each of the giant planets. However- these form 
only a self-consistent set in agreement with the observations and do not 
necessarily 'exclude other models involving a different set of assumptions . 

Figure 10 displays three sets of data in the 7 to 14 ym window in 
the Earth's atmosphere for Jupiter, compared with the calculated emergent 
fluxes for Jim iter Models 3 and 4 of Table 4. It immediately follows 
that the abundance of CH^ above the reflecting layer must be 
v 150 - 200 m-amagat. Model 3 (50 m-amagat of CH^) is clearly too cold 
in the 8 - 14 ym region, while Model 4 (250 m-amagat of CH^) is marginally 

acceptable, but slightly on the warm side. Encrenaz ‘et aL. (1976) have 

_1 

measured this region with a resolution of <L0 cm and found the band 
of CH^ to appear as an emission feature with a temperature of vl45° K. 

This value agrees quite well with other observations and supports a 
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mixing ratio of ^2 x 10 which is the value Orton (1975b) pro- 
poses . The lower abundances of CH^ derived, by de Bergh et al. (1976) 
for the 3v^ band of CH^ perhaps should be re-examined especially with 
regard to the placement of the continuum. 

Other model comparisons suggest that the ammonia particle radius 
is v 3 ym or smaller, which agrees with Kawata and Hansen (1976) who found 
the particles to be less than 2 ym. The larger radii model depressed the 
brightness tanperature in the 9 - 12 ym region approximately 2° K, while 
the 1 -ym-sized particles had little effect in this region. 

The thermal, radiation (15 - 40 ym) calculated in Figure 11 also 
supports the 150 m-amagat CH^ model. Model 4 is seen to be slightly warmer 
than the observed values , although not significantly greater than the 
spectrum of Houck et al. (1975) . From this- comparison we find that one 
suitable model atmosphere which seems to be consistent with all of th e data 
we have presented is that displayed in Table 9. As mentioned earlier, there 
could conceivably be a wide range of other comparably acceptable atmospheric 
models which match the data equally well. 

The effect of phase angle was studied briefly (Figure 12) . The 
band of CH^_ shows almost no variation from center to limb, while the band 
of NH 3 shows a definite limb brightening. The broad band, 10 - 20 ym, was 
found to exhibit limb brightening as was observed with Pioneer 10 and 11 
(Orton, 1975a) . 
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TABLE 9 

An Acceptable Jovian Model Atmosphere 


H 2 

v65 km-amagat 

He/I^ 

0.125 

ch 4 

vl50 m-amagat; ^2 x 10 ~ 

mh 3 

vl3 m-amagat ; NH^/CH^ ml/ 3, 
deep ‘in the troposphere 

Ice 

m3 pm radii 


vl27° K 


Sagan and Salpeter (1976) have proposed a variety of molecular 
and aerosol constituents in the jovian atmosphere. Our models, with 
Tg v200° K in the 5 pa region, are quite consistent with a dense reflecting 
layer cloud at v200° K, xdaich is the appropriate level for a NEL^SH cloud. 
They have suggested ethane, with a mixing ratio v5 x 10" . Also 

suggested are complex organic molecules /particles that reside in the lower 
troposphere, below the ifflU cloud. Our Model 7 has attempted a cursory 
examination of these two additional atmospheric constituents. The 
ratio was 5 x 10~^ and the complex organic particles were assumed to reside 
in the layer below the cloud base and have an average radius of 3 pm. 
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The spectral characteristics were approximated by Kharma et al. (1974) 
for and Khare and Sagan (1973) for the Sg/organics. "The anergent 
flux was essentially unchanged except for the band in scission 

with T-g ^120° K without C^S .^ . This would correspond to Tg ^135° K for 
the higher methane abundance model. We were unable to discern any signifi- 
cant difference in the 8 - 14 ym region with the addition of the complex 
organic particles. This is due to two effects: (1) their absorption fea- 

tures are masked by higher level CH^ and NHg, and (2) any scattering is 
also masked by. the overlying NHg cloud. The 8 - 14 ym region does not 
seen to be the optimum region for such observations. 

Saturn Emergent Flux 

In carparing the calculated and observed brightness temperatures 
of Saturn (Fig . 13) , it is evident that there are areas both of good and 
poor agreement. The models all compare favorably with Low and Davidson 
(1969) in the 5 ym region, with Tg ^120° K being essentially the NHg cloud 
top. S im i l a r ly, good agree m ent is found in the 20 ym region. The major 
problen is in the 8 - 14 ym portion of the spectrum, where NHg and NHg-ice 
particles play an important role. 

The spectra for Models 7 and 9 (100 m-amagat of GH A ) fall between 
the spectra of Models 5 and 6 (35 and 300 m-amagat of CH^, respectively) . 

The only .difference the addition of NHg made was the lowering of the bright- 
ness taiperature by several degrees in the 8 - 12 ym region.. However, the 
models are still almost 10% too warm in this region [cf . 93 + jPk for 8-14 
yen (Low, 1964) and 99 + £pK for 10 ym (Mirphy et al. , 1972) ] . 
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In order to natch these observations, there most be a mechanism 
that provides more to the upper reaches of the troposphere. There 
are at least two possible such mechanisms : (1) is slightly super- 

saturated or (2) convection carries the NH^-ice particles to higher alti- 
tudes. A first order estimate of the NH^ required is only about a factor 
of 2 more than is present in Model 9 in agreement with Woodman et al, (1977) , 

The emergent flux in the narrow spectral region of 1280 - 1330 cm "*" 
is virtually independent of phase angle (see Fig. 14). The band centered 
on the methane fundamental at 1306 cm - '*' had a brightness temperature equal 
to the asympotic mesospheric temperature of the model. An observation in 
this region of the spectrum would give a fairly accurate esti m ate of the CH^ 
abundance on Saturn. The low NHL abundance should be seen in a 920 - 970 cm "*" 
limb darkening. The 10 - 20 um broad band is also seen to darken as the limb 
is approached. If there is a mechanism for punping into the upper tropo- 
sphere of Saturn, it should affect the limb darkening in the 8 - 14 um region, 
perhaps even reversing the limb darkening' to brightening. 

The evidence is insufficiently clear to refine our estimates of 
hydrogen and methane abundances beyond the range our models covered. A 
measurement of the brightness temperature in the 1300 cm - '*" region should 
give a reliable asympotic value for the mesospheric teuperature. From this 
value, the abu ndanc e of methane can be deduced. The effective temperature 
of 97° K seems to fit the observations- quite well , with T g = 93° K being 
clearly too cold. Consequently, Saturn, like Jupiter, must have an internal 
heat source that provides approximately 3.5 times as much energy as Saturn 
receives from the sun. 
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Uranus and Neptune Emergent Flux. 

The paucity of observations of Uranus and Neptune rule out 
definitive conclusions at present- This is also true because our 
models indicate substantial methane, clouds on both planets. Our 
understanding of scattering by Choice particles is limited to that 
of non-absorbing Me scatterers of unknown size or distribution. 

The large abundances reported for H2 and CH^ are an- overestimate 
unless our saturation vapor pressure curve for methane has substantial 
errors at the low tanperature extremes found on these planets. Ab- 
sorption features seen in the visible and near infrared, portions of 
tile spectrum must be greatly enhanced by multiple scattering . 

In our models, the average particle radius was 3 ym; this 
size does not seriously effect the long wavelength thermal radiation. 

If the CH^-ice particles are larger, then the shorter wavelength thermal 
radiation will originate higher in the atmosphere , making it cooler , 
the troposphere will warm and the longer wavelength radiation will appear 
at higher temperatures , A comparison of theory and observation for 
Uranus is displayed in Figure 15. An example of an acceptable model 
atmosphere for Uranus and for Neptune is displayed in Table 10. Our 
Model 5 for Uranus matches the observations rather well and it would 
seen fairly easy to detect a rise in the planet’s effective temperature , 
if it is subject to seasonal variations. 

The spectrum of Neptune (see Fig, 16) with an effective tan- 
perature of 48° K clearly appears too cold. The S(0) and S (1) lines 
of the pressure induced hydrogen transitions are formed near the top 
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of the thick. methane cloud. If the amount of methane was not 
forced to be maximized, the cloud top would be lower and the ef- 
fective temperature would be higher. A T g of roughly 52° K would 
give a better match to the observations. However, this would re- 
quire an internal heat source for Neptune , as has been suggested by 
Murphy and Trafton (1974) . An alternative could be larger size 
scatter ers . The 10 ym radius particles were seen to raise the tropo- 
spheric temperature on Jupiter, and the same effect should take place 
on Neptune. This will raise the brightness temperatures of the longer 
wavelength radiation at the expense of the shorter radiation, which 
would improve the agreement between the observed and calculated 
brightness temperatures for wavelengths greater than 20 ym. A 
measurement of Tg at a wavelength < 20 ym is very important for Neptune, 
and would distinguish between either large particle scatterers or an 
internal heat source. 

With the extensive methane cloud deck, it makes little sense 
to talk of simple reflecting layer models , There appears to be only 
about lOX of the observed amounts of H 2 above the cloud deck and the 
GH^/H 2 mixing ratio is governed by the saturation vapor pressure curve 
in the troposphere and by the temperature minimum in the mesosphere. 
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TABLE 10 

Acceptable Uranian and Neptunian Model Atmospheres 



Uranus 

Neptune 

H 2 above CH^. cloud tops 

7 km-amagat 

10' km-amagat 

He/H^ 

0.125 

0.125 

CH^/Lk in mesosphere 

^2 x 10“ 5 

^3 x 10" 7 

T e 

58° K 

"48° K with large 
particles ■ 



v52° K with small 
particles 
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CONCIUSIONS 

The new methane and amnonia laboratory measurements and the 
modeling techniques described here have a range of applications beyond 
those of the present paper, including a more thorough study of the 
brightness temperature in the 5 to 25 pm region and an exploration of 
the consquences of a very large amount of methane on Uranus and Neptune 
as proposed by Danielson (1977) . This exploration has underscored the 
importance of better understanding the complex refractive indices of 
methane and ammonia in the infrared, and of using more realistic band 
formation models in estimating atmospheric abundances , especially for 
Uranus and Neptune, where much of the observed gas opacity must be due 
to multiple scattering in the methane clouds . A range of new observa- 
tions, for example, broad band brightness temperature measurements for 
Neptune short of 20 pm, can play an important role in distinguishing 
among competing models . In particular the results presented here may 
he of use in the analysis of limb darkening (or brightening) and 
measurements of the emergent spectral flux from Jupiter and Saturn (and 
possibly from Uranus) in the Voyager missions . 
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Figure 1. Laboratory spectra of 21 amagats of at 331 cm path length. 
The two spectra shewn are for temperatures of 195° K and 
120° K, at a 0.5 enf^ resolution. The two absorption lines 
in the 900 cm - ^ region are due to trace amounts of in the 
cell. The calculation of the transmission based on Trafton's 
(1967) technique falls within the noise level of the spectra. 

Figure 2. A comparison of previously published 1L, spectra and our calcula- 
tion for the first overtone band. 

Figure 3. Size distribution from Eq. (5) for two values of a and three 
values of b (Hansen and Travis 1974) . 

Figure 4. Temperature profiles for model atmospheres of Jupiter. (The 
numbers on the curves correspond to those of the models in 
Table 4.) 

Figure 5. Temperature profiles for model atmospheres of Jupiter for dif- 
ferent aerosol particle sizes . Model 3 (3 pm mean radius) and 
Model 5 (1 pm mean radius) are essentially identical, while 
Model 6 (10 pm mean radius) is warmer. 

Figure 6. Temperature profile for model atmospheres of Saturn for 

T e = 97° K. (The numbers on the curves correspond to the 
numbers of the models in Table 5.) 


Figure 7. Temperature profiles for model atmospheres of Saturn for 

T = 93° K. 
e 
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Figure 8. 


Figure 9. 


Figure 10. 


Figure 11. 


Figure 12. 


Temperature profiles for Uranus model atmospheres. (The 
numbers on the curves correspond to the numbers of the. models 
in Table 7.) 

Temperature profiles for Neptune model atmospheres . (The 
numbers on the curves correspond to the numbers of the models 
in Table 7 . ) 

Calculated brightness temperatures for model atmospheres of 
Jupiter in the 7 - 14 ym region of the spectrum. Model 3 in- 
vokes small methane abundances , Model 4 much larger ones . The 

models are compared to three sets of observations as shown. 

~rjT 

Resolution is v0.2 ym, or v20 cm . 

Calculated brightness temperatures for model atmospheres of 
Jupiter in the 100 - 900 cm - ^ region. Observations are as 
follows : (a) Low (1965) ; (b) Harper et al . (1972) ; (c) 

Armstrong et al. (1972) ; (d) Low et al. (1973) ; (e) Aitken and 
Jones (1972) ; (f ) Gillett et al. (1969) ; '(g) Chase et al. (1974) ; 
(h) Wright (1976) . Curves are for small (50 m-amgat) and large 
(250 m-a mg at) amounts of CH^, with and without NH^. 

Linb brightening/ darkening of Jupiter for various spectral 
regions based on Rfodel 4 and compared with the results of 
Orton (1975a). 
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Figure 13. 


Figure 14. 


Figure 15. 


Figure 16. 


Calculated brightness temperatures for Saturn are .displayed 
for Model numbers 5, 6, 3 and 10 summarized in Table 5. The 
wide band-pass tsmeratufes are indicated by the horizontal 
lines for Model 5 and dashed lines for Model 6. Additional 
data are from (i) Nolt (1977) ; and (j ) Ward (1977) . 

Limb brightening/ darkening of Saturn for various spectral 
regions based on Model 9. 

Spectra of Uranus from Rieke and Low (1974) [error bar 
with arrow] , with 'data added fr cm Murphy and Trafton (1974) , 

Low (1966) , (error bar at 500 cm'~J and Gillett and Rieke (1977) boss's . 

Also shown are the calculated brightness temperatures for the , 

He, and CH^ models of Uranus: (Model 5) T g = 58° K, (Model 6) 

T = 68° K, simmarized in Table 7 , The wide band-pass calculation 

are indicated by the horizontal dashed lines . 

Spectra of Neptune from Sieke and Low (1974) [boxes] with Murphy 

and Trafton (1974) (error bar at 500 cm ^ and Gillett and Rieke (1977) 

[error bar with arrow] data added. Also shewn are the calculated 

brightness temperature for the Ha, and CH^ model of Neptune; 

T = 48° K. Ihe wide band-pass calculations are indicated by the 
e 

horizontal dashed lines. 
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